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Abstract 

We investigate nonclassical properties of the output of a Bose-Einstein 
condensate in Milburn's model of intrinsic decoherence. It is shown that 
the squeezing property of the atom laser is suppressed due to decoherence. 
Nevertheless, if some very special conditions were satisfied, the squeezing 
properties of atom laser could be robust against the decoherence. 
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I. INTRODUCTION 



Recently, the development of dilute gas Bose-Einstein condensation (BEC) [1,2,3] has 
opened up the study of atom laser [4-7], which is the matter wave analogs of optical 
lasers. The ideal atom laser beam is a single frequency de Broglie wave with well- 
defined intensity and phase. The first realization of a pulsed atom laser was achieved 
with sodium atoms at MIT by coupling a BEC from a magnetically trapped state to 
an untrapped state using a rf pulse [8]. Then, it was repeated with long rf pulses [9] 
and Raman transitions [10]. 

In recent years, much attention has been focused on the problem of nonlinear atomic 
optics. The four matter- wave mixing was realized in the remarkable experiment [11] 
by making use of the optical technique of Bragg diffraction to the condensate. The 
possibility of optical control on the quantum statistics of the output matter wave was 
pointed out in the framework of nonlinear atomic optical optics [12]. Several dynamical 
analysis concerning the nonclassical properties of the output of the trapped condensed 
atoms have been presented in Refs. [13,14,15]. 

On the other hand, there has been increased interest in the problem about deco- 
herence of BEC and atom laser [16-19]. One fact causing decoherence comes from 
the nonlinear interaction between atoms. For a single-mode condensate the significant 
effect of atomic collisions is to cause fluctuations in the energy and hence fluctuations 
in the frequency, thus causing increased phase uncertainty, eventually lead to decoher- 
ence [16]. In Ref. [19], the influence of the decoherence on quantum coherent atomic 
tunneling between two condensates is studied. It is shown that the decoherence leads 
to the decay of the population difference and the suppression of the coherent atomic 
tunneling. In this paper, we would like to investigate the influence of decoherence on 
the nonclassical properties of the output of a Bose-Einstein condensate by adopting 
Milburn's model of intrinsic decoherence [20]. The influence of decoherence on the 
nonclassical properties, such as sub-poisson distribution and quadrature squeezing of 
the atom laser beam is investigated. It is shown that, under very special conditions, the 
atom laser beam may exhibit stationary quadrature squeezing in such a decoherence 
model. 

The paper is organized as follows: In Sec. II we briefly outline the simple model 
describing the output coupling of the trapped dilute condensed atoms, in which the 
nonlinear interaction between the atoms and the quantized motion of atom center 
of mass in the inhomogeneous magnetic field has been ignored. Then, the quantum 
dynamical behavior of this system in the intrinsic decoherence model is discussed by 
making use of Bogoliubov approximation. In Sec. Ill, the influence of decoherence 
on the nonclassical properties, such as sub-poisson distribution, quadrature squeezing 
effect is investigated. In Sec. IV, there are some discussions. 
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II. THE MODEL 

In this section, we consider the output coupling of the trapped dilute condensed atoms. 
For simplicity, the atoms are assumed to have two states, \T) and \F), with the initial 
condensation occurring in the trapped state \T). State |F), which is typically uncon- 
fined by the magnetic trap, is coupled to |T) by a one-mode squeezed optical field 
tuned near resonance with the \T) — > \F) transition. The Hamiltonian of this system 
is given as [14,15] (k = 1) 

H = cj 6 t 6 + w a a t a + fi(a6 t c + a t 6c t ), (1) 

where b(b*) and c(c') are the annihilation (creation) operators of bosonic atoms for the 
untrapped state |F) and the trapped state |T) with transition frequency c^o, respec- 
tively. a(a^) are the annihilation (creation) operators of the optical field with frequency 
uj a . Here Q = ^Ju a /2e V is the coupling constant, and V is the effective mode volume 
and Sq the vacuum permittivity. In the system (1), the atom-atom coupling and the 
nonlinear interaction between the atoms and the quantization motion of atomic center 
of mass in the trapped state by an inhomogeneous magnetic field has been ignored. 
It was shown that the above system leads to an oscillation behavior of the quantum 
statistics between the optical field and the output atomic laser beam [15]. 



In what follows, we outline the basic content of the Milburn model of decoherence. 
Based on an assumption that on sufficiently short time steps the quantum system does 
not evolves continuously under unitary evolution but rather in a stochastic sequence 
of identical unitary transformations, Milburn has derived the equation for the time 
evolution density operator p{t) of the quantum system [20], 

^ = l[eM~ l -H)p(t) exp(itf) - p(t)}, (2) 

where 7 is the mean frequency of the unitary time step. In the limit 7 — > 00, Eq.(2) 
reduces to the ordinary von Neuman equation for the density operator. It is easy to 
obtain the formal solution of Eq. (2) as follows, 

00 

p(f) = E^(*)p(o)4(0, (3) 

fc=0 

where the Kraus operator A k (t) is given by 

A k (t)= { ^e-^eM-* — )- (4) 
Vkl 7 

Obviously, £r=o 4 (t)A k (t) = I. 

We assume that the initial state of system (1) is described by p(0) = |^ r (0))(^ r (0)|. 
Here, |^(0)) = \cx)t ® l < ^ ) (0))s with \cx)t the Glauber coherent state of the operator c 
characterizing the condensed atoms in the trapped state \T); |$(0)) s = |0)f® \0)f 
represents that there is initially no occupying atoms in the untrapped state \F), and 
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the optical field is in the squeezed vacuum state |£) = 5(^)10) . S(g) = exp^a 1 " 2 — £*a 2 ) 

is the squeezed operator. 

Substituting the p(0) into Eq.(3), we obtain 

p(t) = e" 7 * exp(ac t - a*c) ]T ^-M k \0) TT (0\ <g> |$(0)) ss ($(0)|M tfc exp(-ac t + a*c) 

(5) 

where 

M fe = expH-(# + #i)], (6) 

7 

i?o = <^o& b + uo a a)a + f2(aa6^ + a*a'&), (7) 

#1 = n(ac6 t + a t c t 6). (8) 

If |a| ^> 1, we can adopt the Bogoliubov approximation [21], i.e., neglect Hi in Eq.(5). 
Then in the Bogoliubov approximation, the density operator p{t) can be reexpressed 
as follows, 

p(t) = \a) TT (a\ ®p s (t), (9) 

where 

Ps(t) = e-^Y, ^-exp(-^)|0) FF (0| ® |0(e|exp(<^), (10) 

is the reduced density operator describing the subsystem of the untrapped atoms and 
optical field. 

Now, we confine ourselves in the resonance case, i.e., u = uj a = u. We define the 
operators a(k) (at(h)) and b(k) (tf(h)) as 

ikH ikH 
a{K) = exp( )aexp( J, 

7 " 7 

t ikH f ifcif v 
a' (A; J = exp( )a' exp( j, 

7 7 

&(&;) = exp( ^~^° )6exp(— ^~^° ), 

7 7 

6 t (A .) = ex p(^)6t exp(-^). (11) 

7 7 

It is easy to obtain that 

a (jfe) = cos(— )e- <fcw ^a - *sin( — je'^'ft, 
7 7 

a t( fc ) = cos(— )e ifcw ^at + isin(— )e-^- fe ^6 t , 
7 7 

tO' £■(")' 
6 (fc) = C os(— )e~ ik ^b - lS m( — )e-^ +h ^a, 
7 7 
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6 t (fc ) = C0S (^i) e ^/7 fo t + isin (!5li)eW7) fl t ) (12) 
7 7 

where Q' = \a\fl and e _i61 = a/|a|. By making use of Eq.(12), we can express average 
values of arbitrary operator functionals G(a, a*; b, tf) as following 

Tr[G(a, a t; 6, ^(i)] = e -i* £ i2_L T r[G(a(A;), a t(fc); 6(fc), 6+ (/c) ) 1 0> ^ <0 1 ® |£> 

fc=0 

(13) 



III. THE INFLUENCE OF INTRINSIC DECOHER- 
ENCE ON THE NONCLASSICAL PROPERTIES 
OF THE ATOM LASER BEAM 

In this section, we investigate the nonclassical properties of the atom laser in the 
Milburn's model of decoherence. The average numbers and the fluctuation of the 
output atoms as well as the out-state photon can be obtained by making use of the 
Eq.(13) 

N a (t) = Tr(a!ap 8 (t)) 

= {\ + \ exp( 7 te 2 ^' /7 - jt) + \ exp( 7 te~ 2lfy/7 - 7 t)] sinh 2 r, (14) 

N b (t) = Tr(b% 8 (t)) 

= [--\ exp( 7 te 2 ^ --it)-- exp( 7 te- 2 ^'/T - 7 t)] sinh 2 r, (15) 
2 4 4 

AiV 2 (t) = Tr(aWap s (0) _ [ Tr (a^ap s (t))} 2 

= {^+e^[i exp( 7 te 2 ^' /7 )+i exp( 7 t e - 2 ^' /7 )+^ exp( 7 te 4lfy/7 ) + ^ exp( 7 te- 4lfy/7 )]} sinh 2 r cosh 2 r 
4 2 2 8 8 

+{- + — e- 7t [exp( 7 te 4 ^'/7) + e xp( 7 t e - W/7 )] 
8 16 

-— e- 27t [exp(2 7 te 2 ^' /7 ) + exp(2 7 te- 2in ' /7 ) + 2 exp(2 7 i cos(2ft'/ 7 ))]} sinh 4 r 
16 

+{- - — e- 7 *[exp( 7 te w/7 ) +exp( 7 te- w/7 )]}sinh 2 r, (16) 
8 16 

AiV 2 (t) = Tr(bWbp s (t)) ~ [TT^b Ps (t))] 2 

= {- A -^ \\ exp( 7 te 2 ^' /7 )+^ exp^te^ 77 )-^ exp( 7 te W/7 )-^ exp( 7 te- 4 ^' /7 )]} sinh 2 r cosh 2 r 
4 2 2 8 8 

+{\ + ^e^[exp( 7 te 4lfy /7) + exp^e" 4 ^)] 
-^e- 27 *[exp(2 7 te 2in ' /7 ) + exp(2 7 te- 2in ' /7 ) + 2exp(2 7 tcos(2tf/ 7 ))]} sinh 4 r 

+{\ - 4e- 7 *[exp( 7 te 4 ^'/ 7 ) +exp( 7 te- 4 ^'/ 7 )]}sinh 2 r, (17) 
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where r = 2|£|. 

In order to discuss the quantum statistical properties of the atom laser and the 
out-state optical field, we can calculate the Mandel Q parameters defined as [22] 

«,(^« (< = ..») (18) 

Qi(t) < 0, Qi(t) = or Qi(t) > mean the quantum state of the optical field or the 
atom laser field satisfy sub-Poisson, Poisson or super-Poisson distribution, respectively. 
In Fig.l and Fig.2, the Mandel Q parameters Q a and Qb are plotted as a function of the 
time t for two different values of the parameter 7, respectively We can observe that 
both the optical field and the atom laser beam satisfy the super-poisson distribution 
for any t > in the cases with finite values of the parameter 7. 

We now turn to discuss the influence of intrinsic decoherence on quadrature squeez- 
ing of atom laser and the optical field. We introduce four quadrature operators defined 
by 

X[ a) = l -{a + a\ X( ffl) = l(a-at), X? = \{b + &t), X ® = 1(6 - &t) (19) 
These operators satisfy the commutation relation 

[x{°UJ B) ] = i [Xi b \x$»] = i (20) 
which implies the Heisenberg uncertainly relations 

((AXj a) ) 2 )((AXf ) ) 2 ) > 1 ((AX?y)((AX?f) > 1. (21) 

Squeezing is said to exist whenever ((AX-^) 2 } < 1/4, (i — 1, 2), (j = a, b). In order to 
characterize the influence of intrinsic decoherence on the quadrature squeezing of the 
atom laser and optical field, we calculate the following squeezing coefficients [23] 

a) _ ((AX^) 2 )-0.25 

^ ~ 0^5 ' {22} 

where — 1 < < for quadrature squeezing. We express the squeezing coefficients 
as follows 

S[ a) = 2AT a (t)+Re{e- 7 *[exp(7te- 2 ^ /7 ) + i exp^te 2t{n '~ lv)/ ~ / )+^ exp(-ft e - 2t{n ' +uj)/r )] sinhr cosh re^}, 

(23) 

= 2A^ a (t)-Re{e- 7 *[exp(7te~ 2 ^ /7 )+^exp(7te 2l(n '^ )/7 ) + i exp( 7 te- 2l ^' + ^/ 7 )] sinhr cosh re*}, 

(24) 

S[ b) = 2N b (t) + Re{e- 7 *- 2ie [-exp( 7 te- 2 ^ /7 ) 
+ - exp( 7 te 2i ^'^)/ 7 ) + - exp(7te- 2i ^' +w )/ 7 )] sinhr cosh re*}, (25) 



S? = 2N b (t) - Reie-^-^i-exp^te- 2 ^) 
+ I exp( 7 te 2i(n '^ )/7 ) + \ exp( 7 te- 2i(fy+w)/7 )] sinh r cosh re 10 }, (26) 

where = |. It is obvious that + = AN b (t), and N b (t) is always non- 
negative. So, we need only investigate or to explore the squeezing property 
of the atom laser. In what follows, it is assumed = 9 = 0. We start our analysis 
of the squeezing properties of both the atom laser and the optical field in the limit 
case with 7 — > 00 and finite values of uj and Q', which means not any decoherence 
is presented. Then, the Eqs. (23-26) reduces to the results in Ref. [15], in which the 
squeezing coefficients of both the atom laser and the optical field exhibit the Rabi-like 
oscillation. If 7 is a finite value but remains large, i.e. uj/j <C 1 and Cl'/j <C 1, the 
expressions of the squeezing coefficients and Sf } (i = 1,2) can be approximated 
as 

S{ a) « [1 + cos2fi'texp(-2fi /2 t/7)] sinh 2 r + {cos2o;iexp(-2a; 2 t/7) 

+- cos[2(tf-w)t] exp[-2(fi'-cj) 2 t/7]+i cos[2(fi'+cu)t] exp[-2(fi'+cu) 2 t/7]} sinh r coshr, 

(27) 

S^ a) « [1 + cos2fi'texp(-2fi /2 t/7)] sinh 2 r - {cos 2cut exp (-2^/7) 

+ i cos[2(tf-w)i] exp[-2(fi'-cj) 2 t/7]+i cos[2(fi'+cu)t] exp[-2(fi'+cu) 2 t/7]} sinh r coshr, 

(28) 

Sf } « [1 - cos2fi'texp(-2fi' 2 t/7)] sinh 2 r + {- cos2wtexp(-2cj 2 t/7) 

+- cos[2(tf-u)t] exp[-2(fi'-cj) 2 t/7]+i cos[2(fi'+cu)t] exp[-2(fi'+cu) 2 t/7]} sinh r coshr, 

(29) 

« [1 - cos2fi'texp(-2fi /2 t/7)] sinh 2 r - {- cos 2cjtexp( -2^/7) 

+i cos[2(Q'-o;)i] exp[-2(fi'-cu) 2 t/7]+^ cos[2(fi'+cu)t] exp[-2(fi'+w) 2 t/7]} sinh r coshr, 

(30) 

From Eqs. (27-30), we find that all of the squeezing coefficients and sf^ (i = 1,2) 
tend to a fixed positive value sinh 2 r as the time t — > 00, except that they tend to 
sinh 2 r ± I sinh r cosh r in the special case with f2' = uj. With the further decrease of 7, 
the oscillatory behaviors of the squeezing properties of both the atom laser and the op- 
tical field become frozen. In Fig.3, the squeezing coefficient is plotted as a function 
of time t for three different values of parameter 7. With the decreases of parameter 
7, we can observe rapid deterioration of the Rabi-like oscillation of squeezing coeffi- 
cient. In Fig. 3(b) and Fig3(c), if the time t become very large, the squeezing coefficient 
tends to be a fixed positive value 0.093, which means the decoherence eventually 
completely destroy the squeezing property of the atom laser. However, under very 
special conditions that Q' = uj <C 7, < tanhr < | and 7 is large but remain finite, 

will tend to be a negative value as the time t approaches to infinite, which means 
the stationary state of the atom laser gets squeezed. This can be clearly seen from Fig.4. 
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Recently, the sensitivity of quantum systems that are chaotic in a classical limit 
to small perturbations has been investigated in Ref.[24], and the relation between the 
sensitive and decoherence has been discussed. In what follows, we briefly investigate in- 
fluence of small perturbation on the quadrature squeezing of atom laser. From Eq.(26), 

we can observe that the dynamical behavior of is dependent on the values of the 
exponential facts ^m'-^h _ i] ; 7 [ e ±2ifi'/7 _ ^ 7 [ e -2^/ 7 _ 1 ] > and 7 [ e -2»(n'+ w )/ 7 _ !]_ 

In the case with Q! ~ u < 7, the term ^[e^'-^h _ 1] ~ {2i5 - 25 2 /-f) plays a 

dominantly role in the long time dynamical behavior, where 5 = Q' — uo. In Fig. 5, 

we plotted as the function of time t for different small values of 5. Since Q' 

is dependent on the amplitude a of the initial coherent state of condensed atoms in 

the trapped state |T), we expect that the squeezing dynamical behavior of the atom 

laser is sensitive to the initial state of trapped condensed atoms, just as shown in Fig.5. 



IV. DISCUSSION 

In this paper, we investigate nonclassical properties of the output of a Bose- 
Einstein condensate in Milburn's model of intrinsic decoherence by making use of 
Bogoliubov approximation. It is shown that the squeezing properties of the atom laser 
field can be destroyed by the decoherence under most conditions. However, under 
some very special conditions, the squeezing properties of atom laser is robust against 
the decoherence. This phenomenon is highly dependent on the particular modelling of 
decoherence. 

In what follows, we briefly discuss the relevance of our theoretical results to the 
realistic experimental conditions. By making use of the preliminary atom laser experi- 
ments [8], as a guide to realistic parameter values, we choose uj = 300kHz, Q = 60kHz, 
the initial trapped atom number N = \a\ 2 = 10 6 , the mean frequency of the unitary 
time step 7 = 10 5 MHz, and r = 1 as a example to illustrate our theoretical results. In 
such a condition, the Rabi-like oscillation behaviors of the numbers of both the output 
atom and output photon become frozen after 100/xs. However, the atom laser persists 
in exhibiting oscillation of squeezing coefficient and its squeezing property is finally 
destroyed by decoherence till about 0.15s. 

Milburn's model of intrinsic decoherence is currently a very active field of research. 
Nevertheless, so far little attention has been paid to its study in the context of ultracold 
quantum gases. However, due to their high sensitivity to decoherence, such systems 
could provide an interesting testing ground for the Milburn model. In this paper, we 
provide a first step in this direction by studying the influence of intrinsic decoherence on 
the output properties of a simple atom laser model. A thorough discussion of all other 
sources of decoherence that could be present in an experiment and comparison with the 
results presented in this work should be interesting and necessary in the future study. 
It is also very interesting to investigate the entanglement between output atoms and 
output photons and discuss the influence of intrinsic decoherence on entanglement [25- 
28]. Moreover, it is worth to discuss the possible quantum chaotic dynamical behavior 
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of atom laser and the quantum-classical corresponding by fully taking account of the 
atom-atom interaction. 
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Figure Caption 



FIG.l. The Mandel Q parameter Q a of optical field is plotted as a function of time 
t with ft' — 1 and r = 2 for two different values of the parameter 7, (Solid line) 
7 = 00, (Dot line) 7 = 10 2 . 



FIG. 2. The Mandel Q parameter Qb of atom laser field is plotted as a function of time 
t with ft' — 1 and r = 2 for two different values of the parameter 7, (Solid line) 
7 = 00, (Dot line) 7 = 10 2 . 



FIG. 3. The quadrature squeezing coefficient S% of the atom laser field is plotted as 
a function of time t with 00 = 0.1, <fi = 0, 9 = 0, ft' = 71 and r = 0.3 for three 
different values of the parameter 7, (a) 7 = 00, (b) 7 = 10 3 , (c) 7 = 10 2 . 



FIG. 4. The quadrature squeezing coefficient S2 of the atom laser field is plotted as 
a function of time t with u = 10, <p = 0, 6 = 0, fi' = 10, r = 0.3 and 7 = 10 2 . 



FIG. 5. The quadrature squeezing coefficient S2 of the atom laser field is plotted as 
a function of time t with 00 = 10, <fi = 0, 9 = 0, r = 0.4 and 7 = 10 2 for four 
different values of fi': (Solid line) fi' = 10, (Dot line) Q' = 10 + 10~ 7 , (Dash Dot 
line) ft' = 10 + 2 x 10" 7 , (Dash Dot Dot line) ft' = 10 + 3 x 10~ 7 . 
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Figure 2: The Mandel Q parameter Qb of atom laser field is plotted as a function of 

timp t with O' = 1 anrl r = 9 for two rliflfprpnt valnps nf ttip narampt.pr fSnliH linp^ 
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